Electric fields and valence band offsets at strained [111] 

het er o j unct ions 

Qs . S.Picozzi and A. Continenza 

O ■ INFM - Istituto Nazionale di Fisica della Materia 

D ; 
1^ . Dipartimento di Fisica 

m . 

Universitd degli Studi di L'Aquila, 67010 Coppito (L'Aquila), Italy 

^ . and 

Th '. 

A. J. Freeman 

C^ . Department of Physics and Astronomy and Material Research Center 

Northwestern University, Evanston, IL 60208 (U.S.A.) 

O 

Abstract 
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CS . [Ill] ordered common atom strained layer superlattices (in particular the 

(^ . common anion GaSb/InSb system and the common cation InAs/InSb sys- 

0^ ■ tem) are investigated using the ah initio full potential linearized augmented 



C^ 



O 



% 



plane wave (FLAPW) method. We have focused our attention on the po- 



i-rt ■ tential line-up at the two sides of the homopolar isovalent heteroj unctions 

a: 

O ■ considered, and in particular on its dependence on the strain conditions and 



on the strain induced electric fields. We propose a procedure to locate the 
interface plane where the band alignment could be evaluated; furthermore, we 
suggest that the polarization charges, due to piezoelectric effects, are approxi- 
mately confined to a narrow region close to the interface and do not affect the 
potential discontinuity. We find that the interface contribution to the valence 
band offset is substantially unaffected by strain conditions, whereas the total 
band line-up is highly tunable, as a function of the strain conditions. Finally, 



we compare our results with those obtained for [001] hetero junctions, 
73.20.D, 77.84, 73.20.A 
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I. INTRODUCTION 

As is well known, the opportunity of tuning the potential line-up as a function of the 
strain has a relevant technological importance and thus has been the subject of many recent 
experimental and theoretical works. ETQ However, most of them focused on CuAu-like strained 
interfaces {i.e. oriented along the [001] direction). In this work we have performed ab initio 
full-potential augmented plane waves (FLAPW)l3 calculations (within density functional 
theoryu) for GaSb/InSb and InAs/InSb systems, in order to examine the valence band offset 
(VBO) at strained [111] heterojuctions. To our knowledge, only recently have CuPt-like 
(i.e. [Ill] ordered) systems begun to be investigatedB'tO and the way the charge readjustment 
occurs at the interface is still far from being clear, mostly due to the non-trivial electrostatics 
involved. 

Although both the [001] and [111] directions are polar (the [110] is non-polar), there 
are important differences between interfaces oriented along these two crystallographic axes. 
First of all, the particular geometry for the [111] superlattice leads to two inequivalent 
interfaces:liil consider the structure for the (GaSb)3/(InSb)3 superlattice shown in Fig. |l| 
(note that the common cation case has the same structure, except for the atomic species 
involved). According to convention, we have placed the interface planes tt^ (tt^) at the 
center (sides) of the interface InSb (GaSb) axial bond (as shown in Fig. |I|). Observe that 
the interface Sb is coordinated with three Ga atoms and one In atom in one interface (tTc in 
Fig. |l]), and with three In atoms and one Ga atom in the other (vr^ in Fig. |l]). Thus, we could 
have, in principle, two different valence band offsets and charge accumulations (cr) at the 
interface planes n^ and TCg, namely ac = —Cs, because of charge neutrality requirements.EJ 
Incidentally, we note that this is not the case for the [001] ordered systems, where the bond 
directions and the atomic coordinations result in exactly equivalent interfaces at the center 
and at the sides of the unit cell. 

Furthermore, since the semiconductor constituents (with zincblende structure) are piezo- 
electric materials,li3yj ^ uniaxial strain along the [HI] direction leads to a strain induced 



piezoelectric polarization Ps-HjoEj Following the notation of Ref. |T5|, we obtain: 

|P,| =P, = -^3 714 64 (1) 

where 714 is the piezoelectric constant and €4 is the off-diagonal strain tensor element. As a 
function of the in-plane strain (e_L = """''""'"' , with agub and aeq denoting the substrate and 
equilibrium bulk lattice constants, respectively) and within the linear regime,li3 the strain 
tensor component can be expressed as: 

Cii + 2 C12 , . 

64 = —2 e_L (2) 

cii + 2 C12 + 4 C44 ^ ^ 

The piezoelectric effect can generate electric fields so that:t3 

P =P, + £oXE (3) 

where P is the dielectric polarization and x is the dielectric susceptibility. Note that for 
strained layer superlattices (SLs) with [001] growth axis, the strain tensor is diagonal, so 
that Pg vanishes.Ej 

In this paper, we discuss results obtained for common-cation and common-anion strained 
layer superlattices grown along the [111] direction and discuss the effects of the strain induced 
polarization on the properties of the superlattices. In Sect. II we give some calculational 
details; in Sect. Ill and IV we briefly review the electrostatics involved; in Sect. V we discuss 
how charge redistributes at the interface and how the VBO can be evaluated. Finally, in 
Sect. VI we summarize our conclusions. 

II. STRUCTURAL AND TECHNICAL DETAILS 

We have studied the SL systems (GaSb)3/(InSb)3 and (InAs)3/(InSb)3 with unit cells 
containing 12 atoms (as shown in Fig. ^. For this we have considered an AC/BC type 
heterojunction (the case of the common anion AB/CB type heterojunction is exactly equiv- 
alent) grown on different substrates: (i) on a bulk AC semiconductor {i.e. agub = o-aC) where 



a AC denotes the AC semiconductor equilibrium lattice constant); {ii) on a bulk BC semi- 
conductor {i.e. Qsub = (iBc)i {ii'i) on a substrate having an average lattice parameter {i.e. 
dsub = |(ctyic + ciBc))- The structural parameters for the different systems considered were 
determined according to the macroscopic theory of elasticity (MTE) and are equal to those 
for the ultrathin 1x1 [111] ordered SLs of Ref. |T^ With this choice, we are not considering 
any internal strainlia which would produce different bond lengthsli^l for bonds parallel to the 
[111] and the other directions.By Of course, we expect that variations in the bond lengths 
may greatly affect the magnitude of the electric fields; very large deviations may even cause 
a sing change. However, all this will not change the basic electrostatics involved and, as we 
will show, the physical quantities we are interested in, such as the VBO. 

We performed all-electron full-potential linearized augmented plane wave (FLAPW) cal- 
culations within the local density approximation (LDA)eI to density functional theory with 
the exchange and correlation potentials as parametrized by Hedin and Lundqvist.til Tech- 
nical details (such as wave function cutoff, special k-points mesh etc.) are equal to those 
used in a previous study and reported elsewhere. H'tll In order to check the convergence of 
the valence band offset as a function of the cell dimension, we also performed calculations 
for (GaSb)4/(InSb)4 and (InAs)4/(InSb)4 SLs grown on the same substrates. 

III. STRAIN INDUCED ELECTRIC FIELDS 

Since they are important quantities for determining valence band offsets, we report in 
Fig. H (Fig. §) the common atom Sb (In) core level binding energies E^j^ for the GaSb/InSb 
(InAs/InSb) [HI] ordered SLs grown on an average substrate. In the upper part of Fig. |^ 
(a), we also report the core level binding energies for 3x3 [001] ordered InAs/InSb grown on 
an average substrate. 

First of all, we point out that while in the 4x4 structures the E'^'l ^^^ perfectly aligned 
on a straight line in the bulk region, the same behaviour is not so obvious to extrapolate in 
the 3x3 structures, due to the smaller bulk region. However, in what follows, we will assume 



the linear behaviour also for these structures and verify the validity of this assumption 
afterwards, when comparing the VBOs obtained. 

From Figs. ^ and ^, we observe that the structures considered show electric fields with 
opposite signs at the two sides of the interface. The slopes of the linear trends of Figs. |^ and 
^ give a rough estimate of the absolute values of the electric fields (|E| = E = ^1^ )- Our 
calculations give about 3x10^ V/m (12x10^ V/m) in both the GaSb (InAs) and InSb sides 
of the common anion (cation) heteroj unction, with an error bar ~ 2x10^ V/m. Note that 
this large error bar is due to the very small variations (in particular for the common anion 
case) of E'^i^ in the bulk region, which are of the same order of magnitude of our numerical 
uncertainty upon the core level binding energies themselves. 

Furthermore, we observe that the plot of the core level binding energies for [001] strained 
layer SLs leads to a constant trend (implying a constant value of E^^), as the z coordinate 
(perpendicular to the interface) is varied within each semiconductor bulk region; this con- 
firms the absence of electric fields in the [001] ordered structures. 

IV. PIEZOELECTRICITY AND POLARIZATION CHARGES 

As already pointed out,t3 [111] ordered strained SLs show electric fields, due to piezo- 
electricity. Thus, elementary electrostatics of dielectric media leads to the following relation 
between the macroscopic quantities in each semiconductor constituent (i): 

D^*) = So E(*) + P(^) = eo E^*) + eoX^'^E(^) + P^^^ (4) 

Due to symmetry properties, the transverse {i.e. parallel to the interface) component of the 
piezoelectric polarization vector vanishes; thus, in the following, we will only consider the 
longitudinal {i.e. perpendicular to the interface) component of the electrostatic quantities 
involved. 

The net charge accumulation at an ideal abrupt interface is related to the macroscopic 
polarizations in the two semiconductor constituents, through the following relation:li3 
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a = -(p(2) _ p(i)) . n (5) 

Since the perpendicular component Dn is continuous across the interface, Eq. (H) gives: 

£oE(i) + P(^) = £oE(2) + P(2) (6) 

so that Eq.(^) can be rewritten as: 

cr = +£o(E(2^-E«)-n (7) 

In addition, from Eq.(^) and using Eq.(^), we obtain an equivalent expression for a: 

a = [+60 (x('^E(2) - ;^(i)E(i)) + P,(2) - P« ] • n (8) 

Note that, in the case of growth on a substrate equal to one of the constituent materials, 
only one side of the interface is strained and gives rise to a non- vanishing strain polarization 
Pg. On the other hand, in the case of growth on an average substrate, Pg is non-zero and 
will have opposite signs in the two different materials. In fact, the piezoelectric constantsli3 
have the same sign for all the materials involved, while the two constituents are in biaxial 
tension and compression, respectively, leading to in-plane strains e_^ and e]_ , with opposite 
signs. 

In order to test the validity of the above description, we have performed the double 
macroscopic average (p )Mt3 of the self consistent total (electrons + ions) charge density 
(p^ = p'^'' + p''"^^ for some of the cases considered. We have estimated the net interface 
charge accumulation from the following relation: 

crscF = IP {z)dz (9) 

where the integration is performed between the two bulk regions (referring to Fig. |1| for the 
3x3 case, between hi and 62)- 

We must point out that the use of the symbol a for the charge accumulation does 
not necessarily imply that this charge must have a surface-like distribution. In fact, the 
electrostatics used so far, holds for any volume distributed charge density. However, as will 



be discussed later, the charge distribution at the interface can be reasonably approximated 
as a planar charge. 

Using Eq. (^ and Eq. (|^), we can check the consistency of our results. In the 
(GaSb)4/(InSb)4 SL grown on an average substrate, for example, we obtain ascF ~ 2 
X 10"^ C/m^ from Eq. (^. Using Eq. (|^) and the electric fields as evaluated by the core 
level binding energies [E'^^^ ^ 2.6 x 10^ V/m in the GaSb side of the hetero junction and 
E^"^' '^ 3.3 X 10^ V/m in the InSb side of the heteroj unction), we find cr ?^ 5 x 10^^ C/m^. 
The agreement between this value and (Tscf is within our numerical accuracy, which is es- 
timated to be 1.8 ■ 10^^ C/m^ (considering the error bar on the electric fields) and shows 
the consistency of our calculations. 

However, if one tries to estimate P^^-* and P^^-* from Eq. (|l]) and then substitutes the 
results in Eq. (H), a value for a ~ 49 x 10^"^ C/m^ (a :^ 18 x 10^"^ C/m^) is found, using 
the available calculatedta (experimentalE^O) elastic, dielectric and piezoelectric constants 
for each material. These values seem to be at variance with the integrated value (Jscf- 
The disagreement can be attributed to the use of bulk lattice constants for a region whose 
thickness is as small as a few monolayers and to the fact that the theoretical and experimental 
const ant st3't3 used in Eq. (^) are not consistent with the present calculation, since the MTE 
structures considered do not properly take into account internal strain effects. Let us remark 
that the inconsistency between these values of a and the one obtained ab-initio will not affect 
the evaluation of the VBO, as will be shown (see discussion below). 

Furthermore, from the steeper slopes in Fig. ^ with respect to those in Fig. ^ we observe 
that the magnitudes of the electric fields in the common cation systems are stronger than 
those in the common anion ones. Since Eq. (0) leads to a net charge accumulation propor- 
tional to the discontinuity of the electric fields, it is reasonable to expect a systematically 
smaller charge, a, in the GaSb/InSb than in the InAs/InSb systems, as confirmed by our 
self-consistent results obtained from Eq. (0). 



V. VALENCE BAND OFFSET 

We recall that the potential line-up can be estimated following the procedure used in 
XPS experiments and widely adopted in all-electron calculationsjH^lHa which takes core level 
binding energies as reference levels. According to this method, the VBO is obtained as 
follows: 

A E^ = Ab +AEb (10) 

Here A 6 is an "interface" term (which denotes the core level energy difference between 
atoms in the two bulk regions, i.e. those denoted by "bi" and "b2" in Fig. [^) and A Ei, is 
the so called "bulk" term and indicates the binding energy difference - with respect to the 
topmost valence level - between the same core levels in each semiconductor strained as in 
the heterojunction. This last term takes into account the effects of strain on the electronic 
band structure of each bulk material. 

We must now note that the evaluation of the A b term is usually done assuming that the 
bulk regions of the SL are thick enough so that the atoms taken as reference can be, by all 
means, considered as "bulk" atoms. This procedure implicitly requires that the A b term 
must be evaluated in the limit of infinite distances from the interface, where the reference 
atoms are no longer influenced by the charge rearrangement at the interface. 

It is now obvious that this same procedure cannot be used for the [111] oriented hetero- 
junctions. The presence of electric fields, in fact, makes the A b term vary with the distance 
from the interface, so that its evaluation in the limit of "bulk region" {i.e. infinitely far 
away from the interface) would result in an ill-defined A b. Let us observe, by the way, that 
the same conceptual problem is present for procedures that use the electrostatic potential, 

rather than core levels, as reference energies for the evaluation of the VBO.Ii^ This problem 

no 
has been usually overcomeH'c3 by extrapolating graphically the linear behaviour of the po- 
tential (or the core level binding energy) in the two bulk regions to an "interface plane" , vr, 
arbitrarily taken half-way along the interface bond, and evaluating the difference between 
their intercepts. 

9 



In the following, we illustrate a more rigorous procedure which will also be able to give 
a reasonable estimate of the related error bar. 

A. Interface charge and evaluation of the VBO 

Let us first discuss how the charge is modified at the interface. We can split up this 
charge distribution as the sum of two different effects: (i) a charge accumulation a, due 
to the piezoelectric effect combined with the boundary conditions, which is basically a 
monopole-like distribution and (ii) a charge depletion and accumulation at the two sides 
of the heteroj unction, which is the usual dipole charge of the interface and is due to the 
different chemical nature and electronegativity of the constituents. The charge profile can 
thus be schematically represented as a symmetric monopole centered on some ideal interface 
plane tt plus a dipole which is anti-symmetric with respect to tt. 

According to this scheme and following basic electrostatic arguments, we can now say 
that the monopole charge a leads to a continuous potential across the interface plane. The 
limits for z ^ oo and z -^ — oo of this potential must be the same, since the monopole 
contribution is symmetric with respect to the plane vr. On the other hand, the dipole-like 
charge leads, in these same limits, to a potential discontinuity. Thus, we expect that the 
valence band offset {i.e. the potential discontinuity at the interface) is only determined by 
the usual dipole term and is not affected by the presence of the strain induced electric fields. 

Note that the real total charge density macroscopic average profile (p ) usually has a 
very complex shape, often changing its sign. Therefore, it is not possible to clearly establish 
where tt really is. However, as a rough approximation, we can calculate the center of gravity 

rp rp 

of the positive (j)_^{z)) and negative (j)_{z)) charges: 

z, = l^^ (11) 

! P+{z)dz 

-z^ = l^i^ (12) 

J p_{z)dz 
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where p (z) = p^{z) +p_{z) and the integrals are performed between the two bulk regions 
bi and 62 as in Fig. |l|. (^. We may now assume that the interface region is confined between 
the two extremes z+ and z_. 

We are now able to calculate the A 6 term appearing in Eq. ( |T0|) and therefore the VBO. 
A reasonable estimate for the A b term could then be obtained from the intercepts at z^ 
and Z- of the two straight lines which interpolate the core level energies. The difference of 
these two values can be taken as the error bar. We find, for the systems considered, that the 
values of z^ and z_ are, as expected, not symmetric with respect to the half-bond position 
(which is, in any case, within the (z+,z_) interval), due to two different effects: (z) the 
strain, which modifies the bond lengths at the interface and {ii) the distribution and sign 
of the monopole term a. In all the cases considered, z+ and z_ differ by as much as 2 a.u., 
which is to be compared with bond lengths along the z direction that are usually of the 
order of 5 a.u.. This estimate of the charge centers of gravity confines, with some precision, 
the "interface" to a quite localized region, which is even smaller than a bond length. Within 
this error, we may also say that the interface charge has a surface-like distribution. The 
numerical uncertainty on the A b term can thus be related to the extrapolation procedure 
on the two extremes of the interface region, z^ and z_, leading to a numerical uncertainty 
of about 0.03 eV. 

We may now recall, as already pointed out, that in a [111] ordered SL, we have two 
inequivalent interfaces with, in principle, different A 6's. For each of the structures exam- 
ined, we find that the A b values obtained following this procedure at the two inequivalent 
interfaces differ by as much as 0.02 eV, so that we can consider, as a first approximation, 
A6c = A 6^. 

Furthermore, we also obtained equal A b values for the 3x3 and 4x4 SLs (within 0.02 
eV) , thus confirming the validity of the extrapolation procedure; this also suggests that bulk 
conditions between the two interfaces are recovered even in the smaller unit cell and that 
good convergence for the A b term is reached, as a function of the cell dimensions. Note that 
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the value (0.17 eV) obtained in Ref. ^for InAs/Insb [111] grown on an average substrate is 
in good agreement with the one derived from Fig. ^ (a) (0.19 eV). The small discrepancy 
(about 0.02 eV) is due to the electric fields which modify the core level binding energies and 
which were not taken properly into account in Ref. |^. 

As a last comment, we would like to come back to the structural choices considered. 
In order to ascertain the dependence of the VBO on the internal strain (not considered 
within MTE), we calculated the VBO for superlattices having the experimental internal 
strain displacement: we found a change in the band line-up of less than 0.02 eV (which is 
well within our numerical accuracy). 

B. Results and discussion 

In Table | we report our results for the interface contribution A b, the bulk contribution 
A Ef, and the resulting band offset value A E^j. The superscripts (nr) and (r) indicate a 
semirelativistic and fully relativistic {i.e. spin-orbit coupling included with a perturbative 
approach) treatment of the valence levels, respectively. The final results for A E^^^"' include 
a correction for quasi-particle effects that were omitted in our self- consistent calculations. 



Thus, starting from the results of Ref. ^, we can include self-energy corrections (considered 
uninfluenced by different strain conditions) to the DFT-LDA energy levels by adding the 
quantities A'^^iGaSb/InSb) = -0.03 eV and A'^^{InAs/InSb) = +0.08 eV to our calculated 
values . 

In order to examine the effects of the ordering direction on these quantities of interest, 
we plot our Table I results in Figs. || and | for GaSb/InSb and InAs/InSb interfaces, re- 
spectively, and also report the calculated values for [001] systems with equal pseudomorphic 
growth conditions.^! An inspection of Fig. ^ (a) and Fig. |^ (a) shows that, as a function 
of the substrate lattice parameter, the A b term is almost constant in the GaSb/InSb [111] 
systems (as in all the [001] structures^ - see the dashed lines), whereas in the InAs/InSb 
[111] systems the range in which A 6 varies is slightly larger (about 0.06 eV). This indicates 
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that the charge rearrangement at the interface is quite unaffected by strain conditions, while 
it is much more influenced by the growth direction. In fact, we point out that the A b ab- 
solute value is systematically higher (by as much as 0.15 eV) in the [111] structures than 
in those oriented along the [001] axis; furthermore, we observe that the trend relative to 
the small variations of the A b term as a function of the substrate lattice parameter seems 
to be opposite in the two different orientations of the heteroj unctions; in particular, as the 
substrate lattice constant is increased, the A 6 increases (decreases) in the [111] ([001]) case. 
However, the variations of the A b term are of the same magnitude as the error bar and it 
is therefore impossible to define a clear A b trend as a function of strain conditions. 

Let us now discuss the A Ei, contribution to the band line-up. As previously noticed in 
Ref. ^ the MTE predicts different bond lengths (by as much as 2 %) perpendicular to the 
interface plane, whether we consider the [001] or the [111] ordered heterojunctions. This 
structural difference is probably responsible for the different A Eb term in SLs grown on the 
same substrate (and therefore with the same in-plane strain conditions) but having different 
crystallographic orientation, [001] or [111] (see Fig. ^ (b) and Fig. ^ (b)). 

Finally, we compare the valence band offsets for the CuAu and CuPt-like heterojunctions. 
The perfectly decreasing trend of the VBO as agub is increased, already evidenced for the 
[001] heterojunctions, is still valid for those oriented along the [111] axis (see Fig. ^ (c) and 
Fig. ^ (c)). Furthermore, we observe that the opportunity of tuning the band offset as a 
function of the strain conditions is greater in [111] than in [001] heterojunctions; in fact the 
range in which A E^ varies in going from the smallest to the largest a^ is about 0.8 eV for 
GaSb/InSb and 0.9 eV InAs/InSb [111] interfaces. 

VI. CONCLUSIONS 

We have presented results obtained from FLAPW calculations performed for [111] inter- 
faces, in particular for GaSb/lnSb and InAs/lnSb heterojunctions. In order to examine the 
effects of the strain induced electric fields on the charge redistribution (and hence on the 
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potential profile) near the heteroj unction, we have proposed a simple scheme to locate the 
interface region. To a first approximation, we find that piezoelectric effects lead to a planar 
distribution of polarization charges, which doesn't affect the potential discontinuity. We 
have followed an extrapolation procedure to obtain the interface contribution to the VBO 
which is found to be almost independent of the strain conditions, but is strongly affected by 
the ordering direction. The total potential discontinuity is found to vary by as much as 0.8 
eV for the GaSb/InSb and 0.9 eV for the InAs/InSb [111] interfaces, thus confirming strain 
as an additional degree of freedom to obtain ''ad hoc" band offsets. 
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TABLES 
TABLE L Interface term (A b), strained bulk term (A Eb) and valence band offset (A E^) 

for (GaSb)3/(InSb)3 and (InAs)3/(InSb)3 [111] superlattices as a function of the substrate lattice 

parameter neglecting (A E^^"*" and A Ev ) and including {A Ej^ and AE-J ) spin-orbit effects). 

A Ev "" is the final result (A Ev + A*^ ), to be compared with experimental values (to the 

best of our knowledge, not available at present). Energy differences are considered positive if the 

level relative to the InSb layer is higher in energy with respect to the GaSb (InAs) layer in the 

common-anion (common-cation) system. 





Ab 


A i^f ") 


A 4^) 


A Et^ 


AEi^^ 


A Ei^'"'^ 




GaSb-subs. 


+0.28 


+0.26 


+0.28 


+0.54 


+0.56 


+0.53 


(GaSb)3/(InSb)3 


Av. subs. 


+0.29 


-0.25 


-0.13 


+0.04 


+0.16 


0.13 




InSb-subs. 
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FIGURES 
FIG. 1. Structure for the [111] GaSb/InSb 3x3 superlattice. The atoms denoted by "bi", "b2" 

and "i" indicate the atoms in the two bulk and interface regions, respectively. The dotted lines ttc 

and TTg denote the conventional positions of the interface planes at the center and at the sides of 

the unit cell, respectively. 

FIG. 2. Sb Is core levels binding energies {E^^) (empty circles) in the (a) 3x3 and (b) 4x4 
GaSb/InSb systems (grown on an average substrate) as a function of the z coordinate perpendicular 
to the interface. The energy scale is referred to an arbitrary zero. The dashed lines indicate the 
linear interpolation of the bulk core levels energies; their intercepts with the conventional interface 
planes (full diamonds) define the interface contribution (A h) to the VBO. 

FIG. 3. In Is core levels binding energies {E^i) in the (a) 3x3 and (b) 4x4 InAs/InSb systems 
(grown on an average substrate) as a function of the z coordinate for [111] (empty circles) and 
[001] (empty squares) ordered SLs. Labels and energy scale as in Fig. 2. 

FIG. 4. Interface contribution A 6 (panel (a)), bulk contribution (AE'^ ) (panel (b)) and total 
valence band offset (A Ey ) (panel (c)) for GaSb/InSb interfaces as a function of the substrate 
lattice constant. Filled (empty) symbols and solid (dashed) lines refer to [111] ([001]) ordered 
hetero j unctions . 

FIG. 5. Interface contribution (A6) (panel (a)), bulk contribution {/S.E^ ) (panel (b)) and total 
valence band offset (A Ev) (panel (c)) for InAs/InSb interfaces as a function of the substrate 
lattice constant. Filled (empty) symbols and solid (dashed) lines refer to [111] ([001]) ordered 
hetero j unctions . 
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